function betaq0 = initialguessC(K_c)
if K_c==5
    betaq0 = 1.0e+06*[nan] ;
elseif K_c==6
    betaq0 = 1.0e+06*[nan] ;
elseif K_c==7
    betaq0 = 1.0e+05*[0
   0.000100599368555
   0.000292418883397
   0.000653475186001
   0.001882961528136
   0.004228493468374
   0.016364931003417
   0.039369508921406
   0.060552592788750
   7.154184173029955];

%%%%FOR ALTERNATE KNOT VECTOR WITH NO EXTRA KNOTS IN THE UPPER TAIL:
% % % % % paramspat2 =
% % % % % 
% % % % %    1.0e+04 *
% % % % %                    0
% % % % %    0.000934071112739
% % % % %    0.002471950029104
% % % % %    0.014759376959328
% % % % %    0.042456314638903
% % % % %    0.093144747995384
% % % % %    0.193642931098593
% % % % %    0.996977715667495
% % % % %    1.762557004840515
% % % % %    8.510287122127915


    betaq0 = 1.0e+04*[0
   0.000934071112739
   0.002471950029104
   0.014759376959328
   0.042456314638903
   0.093144747995384
   0.193642931098593
   0.996977715667495
   1.762557004840515
   8.510287122127915];
    
elseif K_c==8
    betaq0 = 1.0e+06*[nan] ;
    
elseif K_c==9
    betaq0 = 1.0e+05*[0
   0.000100599368555
   0.000276706197770
   0.000611735323215
   0.001336804419484
   0.004084372109245
   0.009958743278159
   0.018418164212145
   0.028527204016308
   0.038636296602481
   2.665214171514545
   5.288413921288253] ;  %%%%3 extra knots between the upper two quantile knots.

    
elseif K_c==10
% % % % %     %%%%WHEN WE USE 40 SIMS:
% % % % %     betaq0 = 1.0e+05*[0
% % % % %    0.000093407111274
% % % % %    0.000247136322502
% % % % %    0.007742439076459
% % % % %    0.015858279398802
% % % % %    0.032699876685200
% % % % %    0.090906807236833
% % % % %    0.205394779462923
% % % % %    0.368695135290484
% % % % %    0.565495397560715
% % % % %    1.673458170235159
% % % % %    2.412100023507998
% % % % %    2.782319820154421] ;
    
    %%%%WHEN WE USE 160 SIMS:
    betaq0 = 1.0e+05*[0
   0.000093407111274
   0.000247136322502
   0.008375782491373
   0.017178574297029
   0.031074170169591
   0.089368625090663
   0.204028750136257
   0.367896133208551
   0.565379744304386
   1.670366333221296
   2.407028774153816
   2.775599820154421];
% paramspat1 =
% 
%    1.0e+05 *
% 
%                    0
%    0.000093407111274
%    0.000247081287930
%    0.008370507595604
%    0.017166395413700
%    0.031073648168307
%    0.089368239033921
%    0.204027901633959
%    0.367897779324991
%    0.565383443299138
%    1.670375219510142
%    2.407035442380508
%    2.775599820154421

% MULTIMETHOD:
% 1.0e+05 *
% 
%                    0
%    0.000093407111274
%    0.000266169673236
%    0.008501518201241
%    0.026530791055895
%    0.056510718857519
%    0.119003095521712
%    0.328973569746645
%    0.634070966638323
%    1.029169418154747
%    1.906800378106631
%    2.608010628636406
%    2.959184039453123
    
elseif K_c==11
    betaq0 = 1.0e+05*[0
   0.000083647791545
   0.000233536012222
   0.006675400507612
   0.013308100759681
   0.022392522218964
   0.053466679983981
   0.110361686231627
   0.218872683445784
   0.427865155074241
   0.849806752562495
   1.332025149214577
   1.761298984756446
   1.976067912802103] ;  

elseif K_c==12
    betaq0 = 1.0e+05*[0
   0.000078321560565
   0.000841220551402
   0.004072734188452
   0.009634579525838
   0.016405697795982
   0.032396970215069
   0.068754314763214
   0.127429555563603
   0.234490268801897
   0.382048937965967
   0.648804225468348
   1.717923976701438
   2.431906224989757
   2.789022432680879];
    
elseif K_c==13
    betaq0 = 1.0e+05*[0
   0.000039160780283
   0.000118000084958
   0.000427005270987
   0.004326787023007
   0.009648702700014
   0.016127699701779
   0.031690339217306
   0.068204656871316
   0.127154765829716
   0.234730002723273
   0.383780676130650
   0.743051557550231
   1.575216583455501
   2.172080538564928
   2.984657853904477] ;  %%%%Correlation(ThetaE,ThetaL)
end



end

